Evolution of cosmic structures in different environments in the quasispherical 

Szekeres model 
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This paper investigates evolution of cosmic structures in different environments. For this purpose 
the quasispherical Szekeres model is employed. The Szekeres model is an exact solution of the 
Einstein field equations within which it is possible to describe more than one structure. In this 
way investigations of the evolution of the cosmic structures presented here can be freed from such 
assumptions as small value of the density contrast. Also, studying the evolution of two or three 
structures within one framework enables us to follow the interaction between these structures and 
their impact on the evolution. Main findings include a conclusion that small voids surrounded by 
large overdensities evolve slower than large, isolated voids do. On the other hand, large voids enhance 
the evolution of adjacent galaxy superclusters which evolve much faster than isolated superclusters. 
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I. INTRODUCTION 



At the end of 1970s astronomers provided observa- 
tional evidence that galaxies in the Universe are dis- 
tributed inhomogeneously. Galaxy redshift surveys show 
that the galaxies form structures such as voids, clusters 
and filaments. Although this is only the distribution of 
visible matter, there are some strong indications that vis- 
ible matter does trace the distribution of dark matter, so 
real matter distribution is similar. All these structures 
evolved from small initial fluctuations which started to 
grow after the last scattering moment. However, diffrcnt 
structures evolved in various ways, depending on their 
environment and neighborhood. The present-day den- 
sity contrast [Eq. (flU)) ] of overdense regions is larger 
than 1 [l[ and inside voids it decends to -1 Q. Thus, 
these structures must be described by exact solutions of 
the Einstein equations without such assumptions as small 
value of a density contrast. This paper provides the anal- 
ysis of cosmic structures' evolution which is free of such 
assumptions. The evolution of the cosmic structures in 
diffrent environments is investigated by employing the 
quasispherical Szekeres model which is an exact solution 
of the Einstein field equations. 

The structure of this paper is as follows: Sec. |TT] 
presents the Szekeres model; Sec. IIIII presents the evo- 
lution of pairs of voids and superclusters in the quasi- 
spherical Szekeres model; Sec. IIVI presents the role of 
expansion in the process of structure formation; Sec. |V] 
presents the connection between the results obtained in 
the Szekeres model and the real large-scale structure of 
the Universe. 



II. THE SZEKERES MODEL 

For our purpose it is convenient to use a coordinate 
system which is different from that in which Szekeres 
[3| originally found his solution. The metric is of the 
following form [4|: 



ds 2 



c 2 dt 2 



(e-k) 



dr - <F 



{dp 2 + dq 2 ) 
£2 : 



(1) 



where ' = d/dr, $ = $(t,r), s 
is an arbitrary function of r. 
The function E is given by: 

E(r,p,q) = ^(P 2 + V 2 ) 



±1,0 and k — k(r) < e 



-fp-f^ + C, (2) 



where the functions 5 = S(r), P 
C = C(r) satisfy the relation: 



P(r), Q = Q(r), and 



C 
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e = 0,±1, 



(3) 



but are otherwise arbitrary. 

As can be seen from ([]}, only e = +1 allows the model 
to have all three Friedmann limits (hyperbolic, flat, and 
spherical). This is induced by the requirement of the 
Lorentzian signature of the metric (JTJ) . As we are inter- 
ested in the Friedmann limit of our model, i.e. we expect 
it becomes a homogeneous Friedmann model at a large 
distance from the origin, we will focus only on the e = 1 
case. This case is often called the quasispherical Szekeres 
model. 

Applying metric |1} to the Einstein equations, with the 
assumption that the energy momentum tensor describes 
dust, the Einstein equations reduce to the following two: 
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2M(r) 



k(r) + -A$ 2 (t,r), (4) 
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4TT-^p(t,r,p,q) 



M'(r) -3M(r)E'(r,p,q)/E(r,p,q) 
$ 2 (t, r)[$'(t, r) - *(t, r)E'(r,p, q)/E(r,p, q)] 



.(5) 



In a Newtonian limit Mc 2 /G is equal to the mass in- 
side the shell of radial coordinate r. However, it is not 
an integrated rest mass but rather active gravitational 
mass that generates a gravitational field. Although the 
p function in Eq. ([5]) is a function of all coordinates, it 
can be shown that the density can be decomposed into 
two parts: the monopole distribution and the part which 
has a dipole structure @, 0| : 



(t,r) +e dlp (t,r,p, q). 



(6) 



The function k(r) is another arbitrary function defin- 
ing the Szekeres model. By analogy with the Newtonian 
energy conservation equation, Eq. ((4]) shows that the 
function (— k/2) represents the energy per unit mass of 
the particles in the shells of matter at constant r. On the 
other hand, by analogy with the Friedmann equation and 
from the metric ([1]) the function k determines the geom- 
etry of the spatial sections t = const. However, since k 
is a function of the radial coordinate the geometry of the 
space is now position dependent. 

Eq. (0| can be integrated: 
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2M(r) 



■k(r) 



= c[t-t B (r)} 



(7) 



where ts is an arbitrary function of r. This means that 
the Big Bang is not a single event as in the Friedmann 
models but occurs at different times for different dis- 
tances from the origin. 

As can be seen the Szekeres model is specified by 6 
functions. However, by a choice of the coordinates, the 
number of independent functions can be reduced to 5. 

The equations of motion T a ^ — are reduced to the 
continuity equation: 



p + pO = 0, 

where is the scalar of expansion and is equal to 



(8) 



emerged either due to slower expansion of the space (for- 
mation of overdense regions) or due to faster expansion 
(formation of underdense regions). In the Friedmann 
limit R — > ra, where a is the scale factor and — > 3Hq. 

The Szekeres model is known to have no symmetry [8] . 
It is of great flexibility and wide application in cosmology 
Q and in astrophysics 0, S], and still it can be used 
as a model of many astronomical phenomena. In this 
paper it will be employed to study the evolution of cosmic 
structures in different environments. 



A. Density contrast 

To compare the evolution of different models the 
change in their density contrast is going to be consid- 
ered. Two different types of density contrast indicators 
are taken into account. 

The first one is a usual density contrast defined as fol- 
lows: 



P~ Pb 
Pb 



(10) 



where pb is the background density. 

However, the density contrast defined as above is a 
local quantity and is not covariant with the coordinate 
transformations. The spatially invariant density contrast 
can be defined as follows [101 ] : 



S 



IK 



h a/3 dp dp 



dx a dxP 



K 



dV, 



(11) 



where JeR, and K G R\{0}. This family of the density 
contrast indicators can be considered as local or global 
depending on the size of S. Such a quantity not only 
describes the change of density but also the change of 
gradients and the volume of a perturbed region. So this 
density indicator describes the evolution of the whole re- 
gion in a more sophisticated way than the 6. Here only 
the case I = 2, K = 1/2 will be considered. 

All models presented in this paper are calculated nu- 
merically as, unfortunately, the class of models described 
by an analytical solution is not sufficiently enough to de- 
scribe the considered cosmic structures. However, in the 
e = 1 (and in the e = 0) case the p, q coordinates have 
an infinite range. Therefore, it is more convenient to use 
a diffrent coordinate system where coordinates do have 
a finite range. 



f n n\ — q ^(^' r ) 1. Coordinate system 

&(t,r) — &(t,r)&(t,r)/$(t,r) The surface of constat t and r can be represented by 

<&'(£, r) — $(£, r)E'(r, p, q)/E(r, p, q) a stereographic projection Employing the stereographic 

projection we can map the infinite surface of p, q coordi- 
In the expanding Universe is positive so the density nates to a surface of a sphere which has a finite range of 
decreases. The structures which exist in the Universe, coordinates 9, 4>- 
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Scot ( - ) cos(0), 



After the following transformations, 

p-P 

q-Q 
r 

we obtain 



Scot ( - ] sm(<j>), 



(12) 



time-dependent mass dipole superposed on a monopole. 
However, since S, P and Q are position dependent, the 
axis of the dipole also changes in the space. The func- 
tions S, P, and Q describe the position of this dipole, and 
as can be seen from Eq. (fT6|) . S describes the vertical po- 
sition of the dipole component, while P and Q describe 
its horizontal position. 



B. Model set-up 



-j^(dp 2 + dq 2 ) = (d0 2 +sin 2 



(13) 



The metric ([T) after such transfromations becomes 
nondiagonal: 



2 r 



S' 2 cot 2 



+2S" cot - (Q' sin <p + P' cos < 



(P' 2 + Q' 2 )]}dr 2 



$ 2 



2S'cot - (Q'cos< 



$2 



-P' sin <j>)} drdcj) + 2— (Q' sin 4> + P' cos < 
E 



+S' cot - ) drd<9 - $ 2 (d<9 2 + sin 2 



(14) 



where 



1 

E 



1 — cos 9 



(15) 



and 



To specify model 5 functions of the radial coordinate 
need to be known. Let us define the radial coordinate as 
a value of $ at the initial instant to = 0.5 x 10 6 yr after 
the big bang: 



f := $(r, t ). 



(18) 



However, for clarity in futher use, the ~ sign is omitted 
and the new radial coordinate will be referred to as r. 

Two of these functions will be te(r) and M(r). Let 
us assume that fs(r) = 0. The function M{r) describes 
the active gravitational mass inside the t = const, r — 
const sphere. Let us describe the mass function in the 
following form: 



M(r) = M (r) +6M(r) 



(19) 



where Mq is the mass distribution as in the homogeneous 
universe, and SM is a mass correction, which can be ei- 
ther positive or negative. The SM is defined similarly as 
in the spherical symmetric case: 



SM(r)=47r^ I du<S> 2 (u,to)<P'{u,t )5p(u), (20) 



E' _ S' cos 6 + sin 6 (P' cos <j) + Q' sin < 
~E ~ S 



(16) 



As can be seen, if t = const and r — const, the above 
becomes the metric of the 2-dimensional sphere. Hence, 
every t — const and r — const slice of the Szekeres e = 1 
space-time is a sphere. However, as 5 1 , P and Q are now 
functions of r, the spheres are not concentric. For the 
spheres to be concentric, the following conditions must 
hold: 



P' = 0, 
Q' = 0, 
S' = 0. 



(17) 



Such conditions entail spherical symmetry and the met- 
ric (TT4"]) becomes the line element of the Lemaitre- 
Tolman model [TJ E3- Due to this non-concentricity 
of spheres the density distribution has a structure of a 



where 5p(r) is an arbitrary function chosen to specify 
the SM. Although Sp(r) is not the initial function of 
density fluctuations (since an initial density fluctuation 
is a function of all coordinates) it gives some estimation 
on the initial density fluctuation of the monopole density 
component. 

The next three functions are P(r), Q(r), S(r). All func- 
tions defining the model are presented as each case is be- 
ing considered. The numerical algorithm used to solve 
the Szekeres model's equations is presented in detail in 
Ref. H3. 

The chosen background model is the homogeneous 
Friedmann model with the density: 



Pb = 0„ 



= 0.24 x 



3_fff 

IttG' 



where the Hubble constant is H = 74 km s 1 
The cosmological constant, A, corresponds to J7a 
where tl A = (1/3) (c 2 A/H 2 ). 



(21) 

Mpc" 1 . 
= 0.76, 



density 



III. DOUBLE STRUCTURES 

In this section the evolution of double structures, 
namely a void with an adjacent galaxy supercluster, is 
investigated. Although within the Szekeres model more 
than two structures can be described, such investigations 
of less complex cases may also be useful because they en- 
able us to draw some general conclusions without going 
into too much detail which could easily obscure the larger 
picture. Even then in Sec. [V] a more complex model is 
also investigated and it is found that the rules extracted 
on the basis of the investigations of the double structures 
are still valid for such more complex situations. The evo- 
lution of a double structure was also previously investi- 
gated by Bolejko However, the analysis presented 
in this paper is much more detailed and comprehensive. 



A. Models with P' = = S", Q' =£ 

As mentioned above, if P' = = S' = Q' the Szekeres 
model becomes the Lemaitre-Tolman model. Hence, the 
class of models considered in this subsection is the sim- 
plest generalisation of the spherically symmetric models. 

The double structure of a void and adjacent super- 
cluster can be described in the Szekeres model in two 
different ways. The first alternative is when SM < 0, 
the second when SM > 0. Both these possibilities are 
examined here. 



1. Models specification 
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Model 1: 



Sp = -5 x 1(T 3 x exp[-(r/8Kpc) 2 ] 

5=1, 

P = 0, 

Q = -0.6 ln(l + r/Kpc) 

x exp(-0.003Kpc~ 1 x r). (22) 

Model 2: 

Sp = 1.14 x 1CT 3 x exp[-(r/9Kpc) 2 ] 

5=1, 

P = 0, 

Q = -1.45 ln(l + 0.2Kpc _1 x r) 

x exp(-0.003Kpc _1 x r). (23) 

The density distributions of models 1 and 2 are pre- 
sented in Fig. [1] As can be seen the model with SM < 
has the void in the center, and the supercluster is de- 
scribed by the dipole component of the density distribu- 
tion. In model 2 the converse applies. The overdense 
region is at the origin and the void is described by the 
dipole component of the density distribution. 



FIG. 1: The present-day density distribution, p/pt- Upper 
panel presents model 1 (SM < 0) and lower panel presents 
model 2 (SM > 0). 



2. Evolution 

In this section we compare the evolution of the density 
contrast, S(t,r,6,4>), and the 52,1/2(^,2) density indica- 
tor for models 1, 2, with the corresponding models of 
a single void and with the models of a single superclus- 
ter obtained within the Lemaitre-Tolman model. The 
Lemaitre-Tolman model is considered because within 
this model one can describe a single spherically symmet- 
ric structure. Such a comparison can demonstrate how 
evolution of a structure changes if there is another struc- 
ture in its close proximity. 

Fig. [2] presents the evolution of the density contrast 
of the model 1 in comparison with corresponding mod- 
els obtained within the Lemaitre-Tolman model. The 
Lemaitre-Tolman model was specified by assuming the 
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same condition as the ones in the Szekeres model at the 
initial instant. Namely, the Lemaitre-Tolman model was 
specified by =0 and the profile of the density distri- 
bution. The local density contrast, S, is compared at the 
point of the maximal and minimal density value. Up- 
per panel of Fig. [2] presents the evolution of the density 
contrast inside the void. As can be seen the bahaviour 
of the density contrast in both models is similar. This 
due to the regular conditions at the origin; ie. at the 
origin, where $ = and some other functions are also 
equal to zero, such conditions have to be imposed so 
that there would be no singularity at the origin (for a 
detailed description of the regularity conditions at ori- 
gin see (H). These conditions imply that the origin be- 
haves like a Friedmann model and this is the reason why 
the quasispherical Szekeres and Lemaitre-Tolman mod- 
els are of a very alike evolution pattern at the origin. 
The lower panel of Fig. [2] presents the evolution of the 
density contrast at the very center of the overdense re- 
gion of the model 1 and the corresponding Lemaitre- 
Tolman model. The growth of density contrast in the 
Szekeres model is much faster than in the corresponding 
Lemaitre-Tolman model. The results of this compari- 
son indicate that within the perturbed region of mass 
below the background mass (SM < 0) the evolution of 
underdensities does not change but the evolution of the 
overdense regions situated at the edge of the underdense 
regions is much faster than the similar evolution of iso- 
lated structures. 

The evolution of the density contrast of model 2 (SM > 
0) is presented in Fig. [3J the evolution of the density 
contrast at the point of minimal density is depicted in 
the upper panel of Fig. [3J and the evolution at the origin 
is depicted in the lower panel of Fig. [3J Similarly as 
in model 1, the evolution at the origin in the Szekeres 
model and in the Lemaitre-Tolman model are very alike. 
The evolution of the void, however, is slower within the 
Szekeres model than it is in the Lemaitre-Tolman model. 
This implies that single, isolated voids evolve much faster 
than the ones which are in the neighborhood of large 
overdensities where the mass of the perturbed region is 
above the background mass (SM > 0). 

Now, let us compare the evolution of the S2.1/2 density 
indicator: 



The upper panel of Fig. [4] presents the evolution of 
S2, 1/2 for an underdense region. The lower panel of Fig. 0] 
depicts the evolution of 62,1/2 for an overdense region. As 
can be seen, S2.1/2 f° r the quasispherical Szekeres models 
considered are comparable and the growth of £2, 1/2 f° r 
the Lemaitre-Tolman model is much smaller. This is be- 
cause the volumes of the considered regions are diffrcnt. 
In the Szekeres model the volume is larger than the vol- 
ume in the Lemaitre-Tolman model. 

Figs. [1] and[2]present the shape of the structures with- 
out corrections for the shell displacement. For expample 
the void in Fig. Q] and Fig. [5] (upper panels) seems to 
be almost spherical. In fact this void is squeezed in the 
+Y direction and elongated in the — Y direction [Q 1 7^ 0, 
P' = = S' — see the metric <HU) and Eq. ([IB])]. This 
fact also leads in some regions to density gradients larger 
than in the Lemaitre-Tolman model, hence, such a large 
disproportion in 62,1/2 between the Szekeres model and 
the Lemaitre-Tolman model. 

The results presented above indicate that the evolu- 
tion of the Szekeres model is much more complex than 
the evolution of the Lemaitre-Tolman model. The evolu- 
tion not only depends on the value of the density contrast 
but also on the density gradients and the volume of the 
perturbed region. This is the reason why S 2y i/2 of the void 
in model 2 (SM > 0) is higher than in other models al- 
though the density contrast in this model evolved slower 
then in model 1. Similarly, as can be seen by compari- 
son of Figs [5] and [BJ the overdense region in the model 
with SM < is much larger than in other models, and 
as a consequence the 52,1/2 for this model evolves much 
faster than in other models. The ^2,1/2 provide us with 
information about the evolution of the whole perturbed 
region. 

The evolution of the density at single point is described 
with the local density contrast S. As can be seen the 
evolution of the maximal and minimal density contrast 
depends on the value of SM of perturbed region. The 
evolution of the density contrast inside large and iso- 
lated voids is faster than inside small voids which are 
surounded by highly dense regions. On the other hand, 
the evolution of the density contrast in highly dense re- 
gions in close neighborhood of large voids is faster due to 
faster mass flow from the voids. 



$2,1 



/2 



h a P dp dp 



p 1 dx a dx? 



^/-dctgdrd9d(l), (24) 



Similarly as above two different types of £ are going 
to be considered: 



P> Pb 



C, 



B. 



Models with P' / S' , Q' # 



In this section models of non-constant P, Q and S are 
investigated. The evolution of these models is compared 
with the evolution of models which were considered in 
Sec. UTTA! 



p < p b — > E = V. 

Since the value of Sik depends on units, the results 
presented in Fig. [5] were normalized so they are now of 
order of unity. 



1. Models specification 



Model 3: 



6 



o 
o 



CD 
"O 



CO 

CD 
"O 




3.5 
3 

2.5 
2 

1.5 
1 



0.5 - 



6 8 10 12 14 16 
time [My] 









sz / 




LT 







10 



12 



14 



16 



time [10 a y] 



FIG. 2: The evolution of the density contrast inside the void 
(upper panel), and inside the supercluster (lower panel) for 
model 1 (SM < 0). The curve SZ presents the evolution 
within the Szekeres model; curve LT presents the evolution 
within the Lemaitre-Tolman model. 
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FIG. 3: The evolution of the density contrast inside the void 
(upper panel), and inside the supercluster (lower panel) for 
model 2 (SM > 0). The curve SZ presents the evolution 
within the Szekeres model; curve LT presents the evolution 
within the Lemaitre-Tolman model. 



Sp = -5 x 1CT 3 x exp[-(r/8Kpc) 2 ] 
S= -(r/Kpc) - 4 , 



P = 0.55 (r/Kpc 
Q = 0.33 (r/Kpc) 



0.4 



0.4 



(25) 



Model 4: 



1.14 x 10~ d x exp[ 

0.9 



Sp 

S = - (r/Kpc) 

0.55 (r/Kpc) 8 



-(r/9Kpc) 2 



P = < 



Q = 0.33 (r/Kpc) 



0.8 



(26) 



Fig. [5] presents the comparison of the present-day den- 
sity distribution in models 1 and 3 in colour coded dia- 
grams. It presents the vertical cross-sections of the con- 
sidered structures. The upper panel of Fig. [5] presents 



the vertical cross-section through the surface of <j> = tt/2 
and the lower panel presents the cross section through 
the surface of (f> w tt/6. The comprehensive study of the 
vertical and horizontal cross-sections of similar models 
was presented by Bolejko Fig. [S] also presents the 
vertical cross-sections of models 2 and 4. As can be seen, 
both structures appear to be similar but, in comparison 
with model 1 , in model 3 the dipole component is moved 
down and right. Model 4 on the other hand presents 
the structure moved down and right in comparison with 
model 2. 

The next section discusses the evolution of these struc- 
tures. 



2. Evolution 

The evolutions of the density contrast inside the voids 
and superclusters of models 3 and 1 is very similar which 
needn't be suprising as model 3 has the same <5(r) as 
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FIG. 4: Comparison of 52,1/2 for models with SM > 0, SM < 
0, and corresponding LT model of a void in the upper panel 
and supercluster in the lower panel. 



model 1. Also, the evolutions of the corresponding den- 
sity contrasts of model 4 and 2 is similar. The func- 
tions S, P, Q were chosen so that they reproduce the same 
shape of current structures and the same density contrast 
inside them — that is why that the evolution of a local 
density contrast is comparable for models 1 and 3, and 
for models 2 and 4. However, it is not clear whether or 
not the evolution of 5*2.1/2 is comparable too. When the 
functions S, P, Q are not constant, the axis of a density 
dipole changes. Also, the volume of the perturbed region 
as well as the density gradients can be different. So it 
may be interesting to compare the evolution of the whole 
perturbed underdense and overdense regions of models 1 , 
2, 3, and 4. 

Fig. [7] presents the comparison of S2,i/2 evolution of 
model 1-4. Similarly as in Fig. [3J the values of Six 
were normalized so they are now of order of unity. The 
primed letters denote models of S' ^ ^ P' , Q' ^ 0. As 
can be seen the evolution of 62, 1/2 for all these models is 
also comparable. These results imply that the evolution 
in the quasispherical Szekeres model does not depend on 
the position of the dipole component. As long as the 
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FIG. 5: The present-day colour coded density distribution, 
p/p b . Models with SM < 0. Upper panel — P' = S' = 0. 
Lower panel — P'^O^S'. 



shape and density contrast of the analysed models are 
similar, such models evolve in a very similar way. 



IV. THE ROLE OF EXPANSION 

The faster or slower evolution rate of the previously 
presented models is reflected by their current expansion 
rate. As has been shown above, the evolution does not 
depend on a relative position of the dipole component 
(evolution of models 1 and 3 is similar). Thus, let us 
focus on model 1 and model 2 only. 

Fig. [5] presents the ratio, &sz/@o> of the expansion 
parameter in the considered Szekeres models to the ex- 
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FIG. 6: The present-day colour coded density distribution, 
p/p b . Models with SM > 0. Upper panel — P' = S' = 0. 
Lower panel — P' # S' . 



pansion parameter in the homogeneous background. As 
can be seen, model 1, with SM < 0, has a larger am- 
plitude of this ratio, and the evolution of a supercluster 
in this model is much faster than in the corresponding 
Lemaitre-Tolman model. On the other hand, model 2 
(SM > 0) has smaller amplitude of the 65.2/60 ratio 
and within model 2 the evolution of the density contrast 
inside the void was much slower than in the Lemaitre- 
Tolman model. So clearly the rate of the evolution is 
connected with the rate of the expansion. This conclu- 
sion is also supported by the continuity equation [Eq. 
©]■ 

Still, there remains a question whether the conclusions 
presented at the end of Sec. IIII A 21 about the evolution 
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FIG. 7: Comparison of ^2,1/2 for models with SM < (upper 
panel) and with SM < (lower panel). C corresponds to a 
"supercluster" — an overdense region, and V correspondes to 
a "void" — underdense region. Primes denote models with 
P' # ^ S'. 



of a density contrast being dependent on the mass of 
perturbed region, are not limited to the class of models 
considered in this paper. Are these conclusions general? 
How they are relevant to the real large-scale structure of 
the Universe? These are the questions that are addressed 
in this and the next section of this paper. 

The models presented above were defined by choos- 
ing functions ts, M(r), S(r), Q(r), P(r). As can be seen, 
the functions S(r), Q(r), P(r) describe the position of the 
dipole and even with S and P(r) constant, we are still 
able to reconstruct the cosmic structures. Moreover, the 
functions S, P, and Q are chosen in such a way that they 
reproduce the present day cosmic structures. 

The other functions which are chosen to specify the 
model include M(r) and is(r). Is it also possible to 
choose other set of functions, such as k(r), and tg(r), or 
any other functions such as those described in Ref. [15| • 

However, if we choose for example k(r) = and the 
mass distribution as in model 1 or in model 2, we find that 
the initial fluctuations diminish with time. This is consis- 
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tent with the bahaviour of the Lemaitre-Tolman models, 
where t' B is known to describe the decaying modes and 
k' — growing modes [n| ■ 

Therefore, it is not the SM but the k function that 
is significant for the evolution. This was also noticed 
in spherically symmetric models of structure formation. 
In Ref. [l7| it was concluded that the evolution of cos- 
mic voids is generally driven by the velocity fluctuations 
rather than by the density fluctuations. In Ref. [l8|, 
where the function k could also change in time, the con- 
clusion is similar, i.e. it is this function that plays an 
essential role in the process of the structure formation. 

However, we cannot take any arbitrary k(r) and M(r) 
(or any other function, instead of SM, defining the 
model), because such arbitrary pairs of functions are in 
many cases "unnatural" and lead either to a large am- 
plitude of ts or to a shell crossing singularity (one of 
the conditions to avoid the shell crossing is t' B < — see 
Ref. [4|]), and in most cases — to both of these situations. 
The large amplitude of is undesirable. The function 
ts(r) describes the moment of initial singularity. The ob- 
servations of the cosmic microwave background radiation 
(CMB) indicate that the Universe was very homogeneous 
at the last scattering moment and as a consequence the 
bang time function, ts(r), cannot have a lager amplitude 
than a few thousand years. Larger values of tsif) would 
induce temperature fluctuations on the CMB sky larger 
than observed. On the other hand, models with ts — 
const are known to describe growing modes only [161 ] , so 
the assumption that ts — const seems very natural. If 
we set ts and specify M(r), then k(r) is already spec- 
ified by the Einstein equations. On the contrary, if we 
set ts and specify k{r) then the M(r) is already fixed. 
Therefore, for the class of models which evolve from small 
initial fluctuations and do not have shell crossings during 
the evolution, as well as reproduce structures similar to 
the observed cosmic structures, the conclusions drawn at 
the end of Section UlI A 21 are valid. However, it can now 
be seen that it is not the mass fluctuations that matter 
but the expansion rate. Higher mass in the perturbed re- 
gion slows down the expansion rate — this is a condition 
hindering the evolution of cosmic voids. On the other 
hand, if the mass of perturbed region is below the back- 
ground mass, such region expands much faster than the 
background, leading to the formation of the large under- 
dense regions. Such large voids enhance the formation 
of large elongated overdensities formed at the edges of 
voids, which are usually called walls. 



V. CONNECTION TO THE LARGE-SCALE 
STRUCTURE OF THE UNIVERSE 

The models presented above are models of two struc- 
tures embedded in the homogenueous universe. Although 
the Universe is much more complicated than that, such 
simple models enable us to come to some general conclu- 
sions. In this section the triple structure is considered 
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FIG. 8: The Qsz/Qo ratio. Upper panel presents the ratio 
of model 1, lower panel presents the ratio of model 2. 

and it will be seen that the behaviour of the evolution of 
cosmic structures in this model is similar to this observed 
in the previous models. 

Model 5 is specified by the following set of functions: 

t B = o, 

5 = 1 x 1(T 3 x exp[-(r/20Kpc) 2 ] - 6.5 x 1CT 4 
x exp{-[(r - 35Kpc)/10Kpc) 2 ]}, 
5=1, 
P = 0, 

Q = 0.33 (r/Kpc) ' 8 (27) 

Fig. [S] presents the density distribution. There is an 
overdense region at the origin, followed by a small void 
which spreads to a larger r. At a larger distace from the 
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X [Mpc] FIG. 10: The evolution of a density profile of model 5. The 

profiles correspond to the profile of X — in Fig. [9] Letters 
correspond to different time instants: a — 0.5 x 10 6 after the 

FIG. 9: The present-day colour coded density distribution, big b^g. b 1 5 x 10 9 yr c 5 x 10 9 yr d 10 x 10 9 yr e 

p/p b , of model 5. _ present instant. 



origin, the void is huge and its larger side is adjacent to 
an overdense region. 

This is diametrically different to what happened close 
to the origin where the void adjourns to the supercluster 
only with a narrow cusp. 

The evolution of model 5 is presented in Fig. [TDJ For 
clarity Fig. \TU\ presents only the profile which is rep- 
resented by the X = line in Fig. [51 This profile is 
shown for 5 diffrent time instants. As can be seen, at a 
larger distance from the origin, where the void is large, it 
evolves much faster and exceeds the speed of the evolu- 
tion of the underdense region close to the origin. Another 
significant fact is that the overdense region connected by 
the void across a larger area evolves much faster than the 
supercluster at the origin which is more compact. This 
model exhibits the features of the models previously con- 
sidered. Thus, it might be speculated that the evolution 
of the real structures follows similar patterns. Namely, 
small voids in the Universe which are surrounded by large 
high density regions evolve much slower than the large 
isolated voids. From the perspective of the continuity 
equation the expansion of the space in this region is very 
slow and this is the reason why the voids do not evolve as 
fast as they could. On the other hand, the expansion is 
much faster inside large voids, where the mass of the per- 
turbed region is below the background mass (SM < 0). 
In these situations matter flows from central parts of the 
voids towards the highly dense regions which form at the 
voids' larger sides and enhance their evolutions. In these 
situations matter flows from central parts of the voids to- 
wards the highly dense regions which form at the voids' 
larger sides and speed up their evolutions. 



VI. CONCLUSIONS 

The galaxy redshift surveys show that the Universe is 
patchy with various structures. These structures include 
small voids among compact clusters, superclusters and 
large voids surrounded by large walls or long filaments. 

The evolution of these cosmic structures in different 
environments in the quasispherical Szekeres model was 
investigated. The Szekeres model is one of the most com- 
plex and spatially inhomogeneous exact solutions of the 
Einstein field equations and it has a potential to be more 
widely used in cosmology. Since it is an exact solution of 
Einstein's equations, it enables to investigate the evolu- 
tion of cosmic structures without such approximations as 
linearity and small value of density contrast. Moreover, 
the Szekeres model is flexible enough to describe more 
than one structure. 

Having investigated various models with two or three 
structures within one frame it may be concluded that the 
evolution of the cosmic structures depends on the envi- 
ronment. In perturbed region which mass is below the 
background mass the amplitude of the expansion's fluc- 
tuations is large and as can be seen from the continuity 
equation [Eq. ©], such conditions enhance the evolution 
of cosmic structures. 

The analyses presented in this paper indicate that 
small voids among large overdense regions do not evolve 
as fast as the large voids do. This is because the expan- 
sion of the space is faster inside large voids than inside 
smaller voids. Moreover, this higher expansion rate in- 
side the large voids leads to the formation of large and 
elongated structures such as walls and filaments which 
emerge at the edges of these large voids. 
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